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SUMMARY 

/ Experimental and analytical investigations focusing on secondary atomization 
and ignition characteristics of aluminum/liquid hydrocarbon slurry propellants were 
conducted. Experimental efforts included the application of a laser-based, two-color, 
j forward-scatter technique to simultaneously measure, free-flying slurry droplet 

diameters and velocities for droplet diameters in the range of 10-200 pm. A multi- 
f diffusion flame burner was used to create a high-temperature environment into which a 
dilute stream of slurry droplets could be introduced. Narrowband measurements of 
radiant emission were used to determine if ignition of the aluminum in the slurry droplet 
had occurred. Models of slurry droplet shell formation were applied to aluminum/liquid 
hydrocarbon propellants and used to ascertain the effects of solids loading and ultimate 
particle size on the minimum droplet diameter that will permit secondary atomization. 

For a 60 weight-percent Al slurry, the limiting critical diameter was predicted to be 34.7 

pm, which is somewhat greater than the 20-25 pm limiting diameters determined in the 
experiments. A previously developed model of aluminum ignition in a slurry droplet was 
applied to the present experiments and found to predict ignition times in reasonable 
agreement with experimental measurements. A model was also developed that predicts 
the mechanical stress in the droplet shell and a parametric study was conducted. A 
one-dimensional model of a slurry-fueled rocket combustion chamber was developed. 
This model includes the processes of liquid hydrocarbon burnout, secondary 
atomization, aluminum ignition, and aluminum combustion. Also included is a model for 
radiant heat transfer from the hot aluminum oxide particles to the chamber walls. 
Exercising this model shows that only a modest amount of secondary atomization is 
required to reduce residence times for aluminum burnout, and thereby, maintain 
. relatively short chamber lengths. The model also predicts radiant heat transfer losses to 
\ the walls to be only approximately 3 percent of the fuel energy supplied. Additional 
' work is required to determine the effects of secondary atomization on two-phase losses 
ip the nozzle. 
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INTRODUCTION 


Slurry or gelled propellants in which a solid constituent is suspended in a 
combustible liquid have been of interest for propulsion applications for several decades. 
Depending upon the application, these propellants are advantageous because of either 
their energy content per unit mass or per unit volume. Many solid chemical elements 
are attractive for formulating slurry or gelled propellants. For example, carbon, 
aluminum, and boron all offer substantial increases in volumetric heat of combustion 
compared to hydrocarbon fuels, while on a gravimetric basis, only boron provides an 
increased heating value [1]. For rocket applications, specific impulse and mass density 
control payload capability for a fixed vehicle configuration. Thus, high mass density 
propellant systems such as O 2 /AI/RP-I become attractive possible alternatives to more 
conventional propellant systems. Recent fuel performance and mission studies [2-6] 
have shown that aluminum-based slurry propellants can offer significant potential 
benefits over non-metallized propellants. For example, use of aluminum with O 2 /H 2 
propellants has the potential for increasing the payload delivered to Mars from low earth 
orbit (LEO) by over 20 percent [3]. The addition of aluminum to nitrogen 
tetraoxide/monomethyl hydrazine (NTO/MMH) propellants may also enable planetary 
missions that would otherwise not be possible with neat NTO/MMH [5]. It has also been 
shown that using aluminum can reduce combustion instabilities in liquid hydrocarbon 
propellants [7-8]. 

Several key technical issues impact the use of AI/RP-1 slurries or gels in 
propulsion systems. Among these are rheological properties and ignition and 
combustion characteristics. In the area of rheology, questions concerning the proper 
formulation to provide a stable, easily pumpable slurry with acceptable spray 
characteristics remain unanswered. Recent studies of propellant formulation [9-11] and 
spray characteristics [12] have made inroads in understanding how to deal with the 
conflicting requirements of long-term storability without separation of the solid and liquid 
phases yet maintaining sufficiently low viscosities so that the slurries are easily pumped 
and sprayed. 

Various laboratory-scale experiments have demonstrated combustion of 
aluminum slurry fuels in both rocket thrust chambers [7,8,13,14] and air-breathing 
combustors [15-17]. In all but one of these studies [17], successful ignition and stable 
combustion were achieved. In some cases, combustion efficiencies were as high as 99 
percent. In order to achieve high combustion and thrust efficiencies, the aluminum must 
ignite and burn in the residence time available, and the two-phase losses associated 
with the condensed-phase products must be minimized. Both of these requirements are 
facilitated by having small slurry droplets, since prior work [18,19] has shown that the 
individual aluminum particles in the slurry droplets can form agglomerates that ultimately 
bum as a single aluminum droplet. Good atomization of slurry fuels, however, can be 
difficult to achieve. Fortunately, the additives used to produce desirable rheological 
properties can lead to secondary atomization of large slurry droplets, as observed in 
studies of single droplets of AI/JP-10 [1 1 ,20,21] and other slurries [22-24]. In particular, 
the surfactants used as wetting and dispersing agents promote the breakup of burning 
slurry droplets into a very large number of smaller droplets. The aluminum in these 
small secondary droplets readily ignites and burns. In addition, the particle sizes of the 
aluminum oxide products are smaller, potentially decreasing two-phase flow losses. 
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The focus of the present investigation is this secondary atomization process. 
This report highlights the major objectives and accomplishments of our research. More 
detailed information can be found in previous reports and publications [25-31]. 


PROGRAM OBJECTIVES 

The overall objective of research is the development of a fundamental 
understanding of the ignition and combustion of aluminum-based slurry (or gel) 
propellant droplets using a combination of experiment and analysis. Primary specific 
technical objectives include the following: ■ - ~ ■ ■ 

1 . Understanding the roles of aluminum particle size, surfactants, and gellants 
(or other additives) in promoting or inhibiting secondary atomization of 
propellant droplets; 

2. The extension of previously developed analytical models to address various 
important phenomena associated with secondary atomization that were 
neglected in the simplest model; 

3. Development of a simple one-dimensional model of a slurry-fueled rocket 
engine that incorporates secondary atomization, carrier fuel burnout, and 
aluminum ignition and combustion. This model will allow the evaluation of 
secondary atomization effects on combustion efficiency and two-phase flow 
losses. 

PROGRAM ACCOMPLISHMENTS 
Experimental Systems 

Apparatus and experimental techniques were developed to study the secondary 
atomization of slurry droplets in the 10-100 p.m diameter size range. Specifically, non- 
intrusive optical methods were implemented to measure droplet size and velocity, and to 
detect aluminum combustion of free-flying droplets in a high-temperature environment. 

Particle Size and Velocity Measurements - Slurry droplet and/or agglomerates 
and product particle size and velocity measurements are accomplished using a single- 
particle sizing system based on a two-color forward scattering technique [32]. A 
schematic of the system is shown in Fig. 1. Particle sizing is performed by correlating 
particle diameter with the intensity of near-forward scattering from a particle passing 
through the center of an Ar-ion beam. Right-angle scattering from a He-Ne laser beam, 
coaxially focused inside the Ar-ion beam is used to trigger data acquisition, thus, 
assuring that the particle passed through the center of the Ar-ion beam. A 50/50 beam 
splitter, incorporated in the right-angle coll ection optics, s eparates the 90° scattered light 
signal into two components. The first component triggers the data acquisition, while the 
second is used for aluminum combustion detection, as described below. Particle size 
versus output signal calibration [27] were performed using various diameter pinholes 


2 






□ «g 

CO - 1 


g < z 
? UJ < 
CD 0. 


FIGURE 1 . Schematic of experimental diagonostics for droplet size and velocity measurements and aluminum 
combustion detection. 













(10-1 00 mm) and were found to agree with Mie theory. These calibrations were also 
verified using monodisperse water droplets, the diameters of which were measured 
using video images [27]. Detailed discussion of the particle size and velocity 
measurement system is provided in Refs. [25-27]. 

Aluminum Combustion Detection - An aluminum combustion detection system, 
consisting of a 395 nm bandpass filter, a photomultiplier tube (PMT), and a 10 kHz low- 
pass electronic filter, is incorporated in the 90° scatter optics to detect spectral line 
emissions from vapor-phase aluminum [33,34]. Since vapor-phase aluminum only 
exists during aluminum combustion, these spectral emissions provide a valid indication 
of aluminum combustion. 

Because broadband emission from molten aluminum, AI2O3, and burning 
hydrocarbons produces a low-level background signal on the combustion detection 
channel, a minimum threshold signal level is used to discriminate between valid 
aluminum combustion and this background radiation. From comparisons of combustion 
channel voltages for burning JP-10 droplets and 60 wt% aluminum slurry droplets, a 2 
volt threshold was chosen as the criterion for aluminum combustion since this was the 
highest signal level observed for the JP-10 droplets. Additional details of the system 
development and calibration are presented in Ref. [28]. 

Burner and Atomization Systems - A non-premixed burner produces a 
homogeneous post-flame region that approximates that of a premixed flat-flame burner 
with none of the inherent flashback problems [35, 36]. A schematic of the burner 
system developed is shown in Fig. 2. Slurry is fed to a gas-type atomizer (Spraying 
Systems 1/8J) by a syringe pump. Slurry droplets generated in a spray chamber at the 
bottom of the burner pass through a tube located along the burner centerline, and 
emerge at the center of the burner face. Adjusting the relative composition of the 
N2/CH4 atomizing gas mixture permits matching of the central tube and flame 
stoichiometries. Approximately three percent of the atomizing flow is carried through 
the central tube, while the remainder of the flow is vented. A valve on the vent line 
controls the pressure differential between the chamber and the surroundings, permitting 
a matching of gas exit velocities between the central tube and the burner flames. More 
detailed information concerning the burner and spraying systems and experience in 
dealing with very high viscosity slurries is found in Refs. [25-27]. 

Slurry Droplet Disruption Models 

Analytic models of the processes leading up to droplet disruption were developed 
and applied to aluminum-based slurries. These models included shell formation and 
internal pressure build-up, of which the latter is a major advance in understanding 
secondary atomization of slurries. The models for rigid shell formation, slurry droplet 
ignition, and aluminum burnout were exercised and compared with experimental results 

for droplets with initial sizes ranging from 10-80 pm. A major issue addressed in this 
part of the research is the role of size of the aluminum particles that constitute the slurry. 
Since the size of the aluminum powder used in slurry fuels is of the order of several 
microns [10,1 1], and since the formation of a dense surface layer of particles is one 
mechanism for droplet disruption [21-24,37], the question arises of whether or not 
secondary atomization occurs when the parent droplet contains relatively few aluminum 
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FIGURE 2. Schematic diagram of burner and spray generator. 




particles, rather than thousands. In previous work [20], we studied relatively large 
droplets (200-1200 pm), while in the present investigation, our focus is the secondary 
atomization of much smaller aluminum slurry droplets (20-100 pm). 

Overview - Recent studies [22,24,37] present a simple theoretical framework for 
shell formation, an event which is considered to be an essential precursor to secondary 
atomization. This framework is schematically illustrated in Fig. 3. When a slurry droplet 
is exposed to a hot ambient environment, heat transfer from the gas flow causes liquid 
to vaporize and the droplet surface to regress. As liquid vaporizes at the droplet 
surface, the aluminum particles that were suspended in the vaporized liquid remain 
behind at the droplet surface, causing the number density of aluminum particles at the 
droplet surface to increase. 

If the initial droplet is large enough and there is a sufficient number of aluminum 
particles in the droplet, the particle number density will increase until a rigid shell is 
formed by particles coming into contact with each other. Experimental evidence 
suggests that rigidity occurs when the particle surface layer is nominally three particle 
diameters thick with a solids volume fraction approximately equal to a cubic packing 
(0 S p = 0.524) [22,24], With the rigid shell, further liquid vaporization occurs at a 
constant external diameter. Additional growth of the shell thickness occurs, and mass 
conservation requires that a vapor-void form somewhere in the droplet interior (Fig. 3c). 
Some time interval after rigid shell formation, surfactant pyrolysis makes the shell 
impermeable and liquid vaporization ceases (Fig. 3d) [24,37]. Heat transfer to the 
droplet continues, causing droplet internal pressure and rigid shell stresses to rise until 
the shell fails. 

Limiting Drop Size for Shell Formation - Rigid shell formation requires that there 
are enough solid particles initially present in a slurry drop to form a hollow sphere with a 
specified thickness. The total number of particles in a droplet with an initially uniform 
solids volume fraction is [22] 


Np.tot 


®s,o (^o /d u 



( 1 ) 


The number of particles in the shell is a function of the assumed shell thickness 
and the droplet diameter when the shell becomes rigid and is expressed as [22] 




(d s -2nd u ) 3 ]/dy . 


(2) 


The remainder of the particles are in the liquid slurry core, which has the same solids 
volume fraction as the initial state, i.e., [22] 


N p ,tol - Np, s = 6 S0 Jt(d s -2nd u ) 3 / 6 . (3) 

The size of the droplet after enough liquid has evaporated to form the rigid shell can be 
obtained by substituting Eqs. (1) and (2) into Eq. (3) and solving for d s . This involves a 
cubic equation that is easily solved numerically using Newton's method. 
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FIGURE 3. Schematic of critical shell formation, shell growth, vapor-void formation, and shell sealing. 




The limiting number of particles required to form a rigid shell nd u thick can be 
obtained by recognizing that the smallest possible shell diameter is just twice the rigid 
shell thickness, 


d s,lim “ 2nd u » W 

and using this result, the minimum number of particles to form a shell can be expressed 

Np.lim = 80 sp n . (5) 

The limiting initial droplet diameter can be obtained by setting N Pi tot [Eq. (1)] equal to 
N p ,iim [Eq. (5)] and solving for d 0 : 

1/3 

d o,lim = 2nd u (0 S p/0 s ,o) ■ (®) 


Equation (6) provides a criteria for determining whether a slurry drop of size d 0 
may undergo secondary atomization due to shell formation, provided the constituent 
particle size and initial solids loading are known, thus, for microexplosions to occur, we 
assume that 

d o > d o,lim • ft) 

To evaluate Eq. (6), we assume that 0 S p is that of cubically packed spherical particles, 
while 0 SO is determined from Y s , p s . and pf, i.e.: 


®s,p — 0-524 , 


0 


Pf y £ 


S ’° Pf Y s + p s (l“Y s ) 


(8a) 

(8b) 


The preceding relationships [Eqs. (1)-(8)] were evaluated over a range of 
parameters including those associated with the slurry employed in the experimental 
portion of this investigation. Figure 4 shows d 0 ,iim. expressed as a multiple of 
constituent particle diameter, d u , as a function of Y s . Curves are shown for assumed 
shell thicknesses of 2, 3 and 4 times d u . Experimental evidence [22,37] suggests that 
n=3 is an appropriate value. As expected, d 0 ,iim decreases as Y s increases [37]. The 

experimental 60 wt% Al slurry (d u = 5 pm) was predicted to have a d 0 ,iim of 34.7 pm. 

The dependence of the rigid shell diameter on initial droplet diameter and 
constituent particle size is shown in Fig. 5. Also indicated in the figure are the total 
number of particles and the fraction of these that go into forming the rigid shell. For 

example, an initial 100 pm droplet with 5 pm particles regresses to a rigid shell diameter 
of about 82.6 pm, a decrease of 17.4 percent. A slightly greater percentage decrease in 


8 




= FIGURE 4. Dimensionless limiting initial diameter for the formation of a 

U rigid shell as a function of aluminum mass fraction.. 

•j= 


9 










diameter (21 .3%) is found for a 40 pm droplet. The limiting diameters and number of 
particles are indicated by the vertical lines at the left end of each curve. It is interesting 
to note that even for the largest aluminum slurry droplet investigated (d 0 = 100 pm), 78 
to 90 percent of the total particles are required to form a shell. In contrast, boron and 
carbon slurries typically have much smaller constitutent particles (d u between 0.1 and 1 

pm), and hence, a much smaller fraction of the total particles is present in the shell. To 
illustrate the differences between aluminum and carbon slurries, calculated shell 
diameter and fraction of solids in the shell for a carbon slurry are shown as dashed lines 
in Fig. 5. The carbon slurry had the same mass loading as the aluminum slurry, but d u 

was 0.3 pm. Here we see that very little diameter regression occurs and that the limiting 

initial diameter becomes quite small (d 0 ,iim = 2.15 pm). Based on the results shown in 
Fig. 5, one would expect carbon slurries to exhibit disruptive burning behavior down to 
much smaller initial droplet sizes, compared to aluminum slurries. 

Disruption Model - A conventional droplet vaporization analysis [38], modified to 
account for the solids [37], is employed to calculate droplet diameter, mass, and volume 
up to the critical shell formation time. During this period, 

= -2rh/(jcpfd|) (9) 


with m given by 


c g(T°o Tb) + fyy +c f (Tb T 0 ) 
hfv+c f (T b -T 0 ) J' 

In this particular approach, it is assumed that the surface temperature of the droplet is at 
the boiling point, which is a reasonable approximation. 


m = 


2rcd s k a 
5-1 In 


Assuming constant properties, Eq. (9) can be integrated to yield 


df(t) 



where 


(11a) 


8k 

K = — 9-*n(B + 1) 

Pf c g 

and 

B = c g 

hfv + Cf(Tb-T 0 ) ‘ 


(11b) 


(11c) 


Since only the liquid evaporates, M is simply determined by the volume change, i.e., 
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( 12 ) 


M(t) = Mo-p,*(^-<§(t))/6. 

The shell formation time can be determined by substituting the value of d s when 
the shell is nd u thick into Eq. (11a). This value of d s is readily calculated knowing e s , p 
and 0 s ,o using Eqs. (1)-(3), as discussed previously. 

In the period between t c and t s , we assume that evaporation continues at a rate 
equal to that at the instant of rigid shell formation, as has been done by other 
investigations [22,24,37]. A more exact treatment would have to deal with the 
complexities of flow through a porous medium that is not necessarily saturated with 
liquid [39,40]. In fact, it is only through a detailed treatment that takes into account local 
shell drying and the subsequent pyrolysis of the high boiling point surfactants that the 
processes leading to shell impermeability can be fully understood. Such a treatment is 
beyond the scope of the present work, and thus, the time interval t s *t c is treated as a 
parameter, which is varied to determine its effect. 

During the period between t c and t s , the rigid shell grows in thickness and a 
vapor void forms. With the evaporation rate given by Eq. (10), d p and dy can be 
calculated by integrating simple mass and volume balances. The resulting expressions 
are [22]: 




(d sc -2nd u f - 3e s.° d s.cK(t s *c) 

' s ' c 1 2(e sp -e Si „) 


dj = 3d«K(t s -t c )/2 


(13) 

(14) 


where the rigid shell is assumed to be wet with the liquid fuel. The vapor volume is 
directly determined from Eq. (14) as 

V v = ndJ/6, (15) 

while the liquid volume is the total volume minus the vapor and solids volumes: 


V, = i(dJ c -e Si0 d|-dJ)/6. (16) 

Equations (15) and (16) provide two of the several initial conditions necessary to 
calculate the internal pressure history after the shell seals. 

To gain some insights into the pressure build-up phenomena with a minimum of 
complexit y , the solids-liquid-vapor system is treated as a single thermal lump, with the 
liquid and vapor existing as a saturated mixture of a single component fuel. 

Furthermore, d s , c is assumed to be fixed; hence, total volume does not change as the 
pressure increases. For purposes of illustration, the properties of JP-10 are used. 
Energy conservation for the lumped system is expressed as 
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dU 

dt 


(17) 


q° = -^- = — (M v u v + M f Uf + M s u s ) 


Introducing the mixture quality, the liquid-vapor mixture internal energy is 

U m = M v u v +Mf u f = M m (u f + xufv). (18) 

Treating the vapor phase as a pseudo-ideal gas with a constant compressibility factor, 
Ufv can be related to hfv, as follows: 

Uf, = hfv-P^-Vf). (19) 

For a rigid impermeable shell, the liquid-vapor system mass and volume are constant, 
hence 


v = Vf + xvfy = constant (20) 

Solving Eq. (20) for x and substituting the result into Eq. (18), together with the 
substitution of Eq. (19), yields 


U m = M m 


CpfT+ 


V-Vf 


zRT 


-+v f L 


. D ( zRT 
hfv "P^— p v, 


where we assume Uf = c P fT. 


( 21 ) 


The enthalpy of vaporization as a function of temperature can be approximated in 
terms of T c and T 0 , at which hfv.o is known [41]: 


hfv ~ hfvp 


(1~T/T C ) 

L(1-T 0 /T C )J 


0.375 


( 22 ) 


Since we are dealing with a wet mixture, the pressure and temperature are not 
independent properties. Thus, we employ the Clausius-Clapeyron equation to relate the 
saturation temperature and pressure, i.e. 


P = Aexp(-B/T). 


(23) 


For JP-10, A has a value of 3.069 1 0 9 Pascals, and B has a value of 4704.2 Kelvin [42]. 
Substituting Eq. (22) and (23) into Eq. (18) and differentiating with respect to time yields 
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where 


and 


= M m (a 1 +a 2 +a 3 )— 

(24a) 

a i - c pf > 


(24b) 

( V — Vf ) 


(24c) 

a 2 - 


a 3 = (v- v f )ABexp(-B/T)/T 2 , 

(24d) 

D “0- 375 hfy 0 

p “ /- ~ \ 0.375 i-r- t-\0.625 * 

(249) 


The liquid specific volume Vf has been assumed constant. 

Substituting Eq. (24) into the energy balance (Eq. 17) for the solids-liquid-vapor 
system and solving for the temperature derivative yields 


£[ = qc 

dt [Mm( a 1 +a 2 + ct 3) + ^s c ps 

where it is assumed that u s = c ps T. 

The integration of Eq. (25) provides the temperature history of the slurry droplet 
after the shell seals, and via the Clausius-Clapeyron relation (Eq. 23), the pressure 
history is determined. The instantaneous convective heat transfer rate is given by 

qo = h it d| c (T„ - T) . (26) 

Consistent with a fixed volume during pressure build-up, we assume that the 
microexplosion event is triggered by the failure of a brittle shell. For a thin shell, the 
internal stress level is readily calculated knowing the pressure difference across the 
shell: 


a 


APr 
25 ' 


(27) 
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Model Results - Figures 6 and 7 illustrate 5, Mf/Mf i0 , and a as functions of time 
after rigid shell formation for initial droplet sizes of 50 and 200 pm, respectively. Thus, 
at time equals zero, 5 is 3d u or 15 pm. For the 50 pm drop-size, very little slurry 

remains after the critical shell is formed, so 5 grows only slightly (~ 6%) as the remaining 
solids are added to the shell as the liquid evaporates. In contrast, the shell of the 200 
pm drop has the potential to increase up to 83% in thickness, if evaporation continues 
without the shell becoming impermeable and sealing the remaining fuel inside the shell. 

The shell stress, a, is shown as a family of three or four lines with each line 
corresponding to a different assumed shell sealing time. For example, in Fig. 6, the 
stress build-up is shown for assumed shell sealing times of 0.01 , 0.2, and 0.5 ms. At 
the instant the shell seals, the mass fraction of the JP-10 becomes fixed and the 

pressure, and consequently the stress, begins to increase. For the 50 pm droplet case, 
we see that the stress increases at essentially the same rate, regardless of the shell 
sealing time, as a consequence of the shell thickness being identical for each of the 

three histories shown. For the 200 pm drop, the stress increases more rapidly for the 
earlier sealing times since the shells are thinner and have not yet reached their 
maximum thickness. We note also that stress levels increase with initial droplet 
diameter. This is a consequence of the larger mean shell radius, since the pressure 
levels for all cases ranged from an initial atmospheric value to 2000 kPa. 

The above results imply that larger droplets are more likely to fragment than the 
smaller droplets because of higher shell stresses; and furthermore, the consequence of 
the larger drops' fragmenting is likely to be much more beneficial since the potential for 
spewing out raw slurry, rather than just producing shell fragments, is greater. In fact, for 
the 50 pm droplet, only shell fragments can be produced if sealing occurs later than 
0.01 ms after the critical rigid shell is formed. 


Limiting Drop Size for Disruption - Theory and Experiment 

An experimental study was conducted to ascertain the minimum slurry drop size 
that would allow disruptive burning and the results were compared with predictions from 
the theory discussed above. 

Size and ignition data for the 60 wt% slurry at x=5, 10, and 15 mm are presented 
in Fig. 8-11. Figures 8 and 9 show respective size distributions for a cold and a hot 
flow, 5 mm above the burner. The most obvious feature of Fig. 8 is the apparent 
bimodal particle size distribution. This distribution may be due to atomizing system 
characteristics or might be an artifact of the system collection bias since a PDF has not 
been applied to the data [32]. 

Before comparisons between experimental and theoretical results can be made, 
the axial location at which droplet combustion begins must be determined. Combustion 
is unlikely to start at the burner surface since the droplets must first pass through the 
cool, oxygen-deficient core of a diffusion flame above the central tube. Combustion is 
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FIGURE 6. 50 wt% aluminum slurry droplet (d 0 = 50 pm) shell stress, liquid fuel 

mass fraction, and shell thickness as functions of time, assuming 
various shell sealing time intervals (d u it = 5 pm). 
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FIGURE 9. 60 wt% aluminum slurry particle size distribution 5 mm above the 

burner with flame present. 
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FIGURE 10. 60 wt% aluminum slurry particle size distribution 10 mm above the 
burner with flame present. 
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FIGURE 1 1. 60 wt% aluminum slurry particle size distribution 15 mm above the 
burner with flame present. 


21 




burner face as evidenced by the height of the surrounding burner flames. Figs. 8 and 
instead likely to begin at the flame tip where heat and oxygen are present. Although the 
height of this flame is not clearly defined, the tip is between 3 and 7 mm above the 
9 (x=5 mm), which show no major differences between the hot and cold flow size 
distributions, further support this premise of delayed droplet combustion activity. 

The data in Fig. 10 (x= 10 mm) show a slight shift toward smaller particles in 
comparison to the x=5 mm position. This is particularly visible in the 5-15 pm range. 
Aluminum combustion is first seen to appear in this same size range, as would be 
expected, since faster liquid burnout results in shorter aluminum agglomerate ignition 
times for the small slurry droplets. 

Figure 1 1 (x=15 mm) shows a sharp decrease in the number of particles in the 
20-40 pm range and a substantial increase in the number of particles in the 10-15 pm 
region. These changes do not result from simple liquid evaporation since calculations 
indicate that agglomerates formed by 20-40 pm droplets should give rise to an increase 
in particle number in the 20-30 pm range rather than the 10-15 pm range. These facts 

lead to the conclusion that secondary atomization of particles in the 20-40 pm range 
occurs between the 10 and 15 mm axial locations. This conclusion is also consistent 
with visual observations of sudden bursts of glowing particles in this same region and 
higher in the flame, with very few below 10 mm. 

Average particle velocities for the data in Figs. 8-1 1 were calculated to be 
approximately 7 m/s, ranging from 2-10 m/s. Assuming that particles travel at this 
average velocity and that droplet life histories begin 5 mm above the burner instead of 
at the burner face, as discussed above, then particles will reach the 10 mm axial 
location in 0.71 ms, and the 15 mm location in 1 .4 ms. It is important to note that these 
calculated times-of-flight are highly sensitive to the location at which droplet combustion 
activity begins. With these assumptions, experimental results correlate reasonably well 
with prior theoretical analyses of aluminum agglomerate ignition times [18,19]. 

As mentioned above, Fig. 10 (x=10 mm) shows some aluminum combustion in 
the 5-10 pm size range. The time for a particle to reach this axial location (t=0.7 ms) is 
comparable to calculated agglomerate ignition times. For example, calculations predict 
that a 14.3 pm slurry droplet, forming a 10 pm agglomerate, would require 0.6 ms for 
aluminum combustion to begin, agreeing favorably with the experimental value of 0.7 
ms. 


Since rigid shell formation is a precursor to secondary atomization, calculated 
rigid shell formation times should match observed disruption times if the shell sealing 
interval is short. Calculations show that a 40 pm slurry droplet should form a rigid shell 
in 0.95 ms. This time interval falls between the particle times-of-flight required to reach 
the 10 mm location (t=0.71 ms), and the 15 mm location (t=1.4 ms). The sharp 
decrease in particle number in the 20-40 pm range seen in Fig. 9 (x=15 mm) occurs in 
this same region. However, it is interesting to note that the lower end of this size range, 
20 pm, is somewhat less than the predicted d 0> iim of 34.7 pm. 
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This discrepancy between experimental and theoretical minimum diameters 
could be due to experimental uncertainty in measuring particle size, primarily associated 
with the unknown slurry droplet index of refraction; however, light scattering in the near- 
forward direction is relatively insensitive to particle index of refraction. Another reason 

for this discrepancy may be that d 0 ,iim of 34.7 pm, was calculated assuming a critical 
shell thickness of 3d u [22]. If this critical thickness is instead assumed to be only 2d u , 
then the predicted d 0 ,iim becomes 24 pm, much nearer the experimentally measured 

limiting diameter of 20-25 pm. A difference in critical shell thickness could be due to 
surfactants or particle morphology affecting particle binding. The critical shell thickness 
of 3d u was determined using a carbon slurry (d u = 0.3 pm) containing unknown 
surfactants and stabilizers. If these additives are different in type or concentration from 
our aluminum slurry, the shell formation characteristics may be affected. Similarly, if the 
ultimate particles in the carbon slurry were more uniform in size or smoother in shape, 
the shell formation processes could again be affected [24]. 


One-Dimensional Combustion Chamber Model Development 

Overview - Recent performance analyses of slurry fueled rockets [2-6] predict 
that Al slurry propellants may increase maximum payload over neat liquid systems. 
However, these studies neglect additional radiation and two-phase flow losses resulting 
from solid combustion products and increases in droplet combustion times due to solid 
agglomerate combustion. Accounting for these three factors is critical in accurately 
predicting the effects of using metallized propellants as rocket fuels. 

The research discussed in previous sections of this report focuses on the 
mechanism of secondary atomization as a possible means of reducing combustion 
times and two-phase flow losses. However, this work is devoted to understanding the 
droplet combustion and secondary atomization processes themselves and does not 
involve the effects of secondary atomization on engine performance. Consequently, a 
one-dimensional engine model was developed to provide a peliminary evaluation of the 
effects of secondary atomization, two-phase flow losses, and radiation heat transfer on 
engine performance. Coding a new model for the combustion chamber provides some 
benefits over using a standard code such as TDK. First, the new model will permit the 
inclusion of secondary atomization effects, and second, the model will allow a clearer, if 
simpler, understanding of the physics involved in the problem. Therefore, a numerical 
code incorporating the following effects has been written to model a rocket combustion 
chamber: 

• Multiple droplet size classes 

• Gas-phase chemical equilibrium 

• Two-phase flow losses 

• Hydrocarbon evaporation and burnout 
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• Al combustion with surface condensation of AI2O3 

• Simple secondary atomization 

• Radiation from condensed products to the chamber walls 

Chamber dimentions and propellant flow rates have been taken from Galecki's work 
[43], to provide realistic model inputs. 

Model Description - The propellants used in the engine model are a JP-10/AI 
slurry and a preheated gaseous O2 oxidizer. JP-10, a pure hydrocarbon (C10H16). was 
chosen as the slurry hydrocarbon component instead of RP-1 to avoid the complexity of 
modeling multi-component droplet combustion. 

The combustor flow is modeled using a single product-phase containing both 
gases and small AI2O3 fume particles, and three additional flow phases for each droplet 
size class: a liquid hydrocarbon phase, an aluminum phase, and a phase containing 
large AI2O3 agglomerates. The AI2O3 must be separated into two flow phases because 
of the two oxidation mechanisms inherent in aluminum combustion [44-46]. In the first 
oxidation mechanism, large AI2O3 agglomerates are formed through droplet surface 
condensation/oxidation, and in the second, very small AI2O3 fume particles are 
produced through vapor-phase oxidation. Since the AI2O3 agglomerates are much 
larger than the fume particles and are attached to the surface of the aluminum droplets 
it is necessary to model the AI2O3 agglomerates separately from the fume particles. 
Including the fume particles in the gas-phase flow requires the assumptions of no 
temperature or velocity slip between the fume particles and the gases, but greatly 
simplifies the problem solution. 

As an Al droplet burns, it continually produces AI2O3 fume particles at the oxide 
boiling point which subsequently equilibrate with the gas-phase temperature and 
velocity. This process results in a range of fume particle temperatures and. velocities as 
new particles are created and older particles continue to equilibrate with the gases. If 
the temperature and velocity slip were not neglected, many additional flow phases 
would be required to accurately model this range of fume particle temperatures and 
velocities. A detailed justification of the no-slip conditions is presented in Ref. [29]. 

Mass Conservation - For M droplet size classes, the following equation can be 
written for system mass conservation assuming steady-state conditions and that no 
mass is added to or removed from the chamber except at the injector face and the 
chamber exit: 


drhg = dm LH + drh A | + dril Ai 2 o 3 ^ 
dx j = | dx dx dx ). 


(28) 


In Eq. 28, x is axial position in the combustor; j, a particular droplet size class; and m g , 
the gas mass flux. The liquid hydrocarbon, aluminum, and AI2O3 agglomerate mass 
fluxes for a particular droplet size class are given as rh|_n, m A |, and m A | 2 o 3 . 
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respectively. The terms on the right hand side of Eq. 28 are found from a hydrocarbon 
droplet gasification model and an aluminum droplet combustion model. The 
development of all three terms is similar and is outlined below. 

The change in mass of a single slurry droplet in a size class, j, during a time 
interval dt is: 


m change,j 


dm drop,j dt 
dt 


(29) 


where m C hange,j represents the hydrocarbon, aluminum, or AI2O3 mass change. The 
value of dmdrop/dt for the liquid hydrocarbon is found from the hydrocarbon gasification 
model while dm’drop/dt for the aluminum and AI2O3 are both determined from the 
aluminum combustion model. Relating dt to dx through the chain rule and the velocity 
relationship 


j p,i 


dx 

dt 


(30) 


yields 


m change,j 


dm drop,j 

dt 



(31) 


where u P j is the velocity of droplets in the j th size class. This velocity is calculated 
through a drag analysis developed later in the report. 


Equation 31 can be converted to a change in mass flux by multiplying m C hange,j 
by the total number of droplets in the j th size class passing through the control volume 
per unit time, N/xj. Substituting this result, Eq. 28 becomes 


dm g 

M 

= -I 

N 

1 

dm drop "1 

+ 

f dm drop ^ 

+ 

dm drop > 


► 

dx 

j=i 

X 

j u P.j 

dt J 

LH 

l dt J 

Al 

^ dt j 

ai 2 o 3 

j. 


The value of N/xj for each droplet size class can be determined from the total 
initial slurry mass flux and normalized droplet size distribution. Given an initial 
normalized droplet size distribution, it is necessary to determine N/xt 0 tal before N/xj can 

be found for the individual size classes. N/xtotal can be found from the following mass 
conservation expression: 


nr 'slurry,o ~ 


N 


total j=f 6 ' ) 


(33) 
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where Nj represents the percentage of total droplets in a given size class, given by the 

normalized size distribution. The variables, dj t o and po, represent the initial diameter of 
the j th size class and the initial slurry density, respectively. With N/x tota | known, N/xj 
for each size class can be calculated from 


N 

x 


lj 


Nj — 

1 X 


It ota I 


(34) 


Neglecting droplet breakup and secondary atomization, N/xj is constant throughout the 
chamber. This leaves the following three variables for each size class to be determined 

from the droplet gasification models: (dm drop j /dt )|-H' ( c,m drop.j / dt ) A |’ 

(dm droPi j / dt ) Al203 • 

Fnergv Conservation - The steady-state energy balance for the system, including 
radiation heat losses, can be expressed as 


d(mh)g 

dx 


M d(mh) LH 
ft dx 


d(rhh) A[ 

dx 


d(mh) A | 2 o 3 


dx 


-27tRq r , 


Jj 


(35) 


where h represents total specific enthalpy (h C hem + hsens). R *s the chamber radius, and 
q r is the radiation heat flux from the flow to the wall. Since gas flow optical properties 
are dominated by the small AI2O3 fume particles [47], the radiation term is independent 
of the individual droplet size classes and is not included in the summation term. 


Each of the bracketed terms on the right hand side of Eq. 35 can be expanded as 


d(rhh) 


dx 


. dh .dm 
m— +h— 
dx dx 


J j 


(36) 


Substituting for dm/dx yields 


d(mh) 


. dh , 

m _ X h 

N 1 d^drop 


dx 

j 

111 — ■ — T 1 1 

dx 

l TU P dt j 

-M 


(37) 


Because the enthalpy required to heat the liquid hydrocarbon from 300K to the 
hydrocarbon boiling temperature (438.9 kJ/kg at P=6 atm) is comparable to the 
hydrocarbon latent heat of vaporization (-286.7 kJ/kg at P=6 atm), hydrocarbon heat up 
must be accounted for in the hydrocarbon vaporization model presented later in this 
report. In the vaporization model, the bulk droplet temperature is assumed to remain at 
the initial droplet temperature. Heat transfer from the gas heats a thin surface layer of 
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liquid hydrocarbon from the initial temperature to the boiling point and vaporizes this 
hydrocarbon layer. Using the constant bulk droplet temperature assumption yields a 
constant liquid hydrocarbon specific enthalpy, which reduces Eq. 37 for the liquid 
hydrocarbon to 




IS 


d(mh)|_H 


h LH,T 0 


N 1 dm drop LH 


X Ur 



( 38 ) 


where h|_H,T 0 is the liquid hydrocarbon specific enthalpy at the initial droplet 
temperature. 

Following hydrocarbon burnout, an agglomerate of aluminum particles remains 
[18,19]. The aluminum agglomerate temperature rises from the hydrocarbon boiling 
temperature, through the aluminum melting point, to the aluminum boiling temperature 
as heat is transferred from the gas flow to the agglomerate. The enthalpy required for 
this temperature increase is significant compared to the total system enthalpy and the 
aluminum enthalpy of vaporization, hfv,Ai- Consequently, agglomerate ignition/heat up 
should be modeled. The agglomerate ignition/heat-up process is approximated as 
convective heat transfer from the gas flow to a spherical, uniform-temperature 
agglomerate. Heat transfer to the agglomerate causes the agglomerate temperature to 
rise from the hydrocarbon boiling temperature to the aluminum melting temperature. 
The agglomerate temperature is then held constant until sufficient energy to melt the 
entire agglomerate has been transferred from the gas-flow, at which point aluminum 
combustion begins. For the heat-up and phase-change processes, the general energy 
conservation equation (cf. Eq. 37) using the above assumptions reduces to 


d(mh) A( | [ . dhl 
dx || " L m dxj) 


J_r A p h(T g -Tai) ^ 

X Up ^ m drop A| c p,AI J j 


(39) 


where h is the average heat-transfer coefficient. 

After the agglomerate phase-change, Eq. 37 for the aluminum mass flux 
becomes 







d(mh) A | 




N 1 dmdropAi 
x u D dt 


(40) 


where h A |j m A| is the liquid aluminum specific enthalpy at the melting temperature. 

Since AI 2 O 3 agglomerate forms on the aluminum droplet surface, the AI 2 O 3 
agglomerate is maintained at the aluminum boiling temperature as long as any 
aluminum remains in the droplet. For the time period prior to aluminum burnout, Eq. 37 
reduces to 
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d(mh) A | 2 p 3 

dx 


h AI 2 0 3 ,T b ,AI 


N 1 dmdrop.AlpO 


£2. 


X Ur 


dt 




for the AI2O3 agglomerate, where h A i 2 o 3 ,T b A) is the AI2O3 specific enthalpy at the 
aluminum boiling temperature. After aluminum burnout, the AI2O3 agglomerate mass 
remains constant but heat transfer to or from the gas flow causes the agglomerate 
temperature to change. Following the development for aluminum heat up (Eq. 39), Eq. 
37 reduces to 


d(mh) A | 2 o 3 


r . dhl 

lYI 


. N 1 

m 

Aph(T g -T A | 2 o 3 ) ] 

dx 

i 

L dxj 

j 

1 1 1 

X Up 

^ m drop,AI 2 0 3 c p,AI 2 0 3 > 


(42) 


for the AI2O3 agglomerates during the time period after aluminum burnout. The total 
system energy balance (Eq. 35) is represented by appropriate substitutions of Eqns. 38- 
42. 


Radiation Heat Transfer - Radiation from the solid combustion products to the 
chamber wall is a participating medium phenomena requiring the solution of the 
radiative-transfer equation. Given a scattering albedo, the transfer equation can be 
solved numerically, but requires a great deal of computer time. Fortunately, the transfer 
equation can be simplified and solved using one of several approximations based on the 
value of the optical thickness, Kd, which is define as 

s s 

k D = J(a + a s )ds = Jicds , (43) 

0 0 

where a and a s are the absorption and scattering coefficeints, respectively, and S is the 
characteristic radiation path length. If k, defined as (a+o s ), is constant along this path, 
Kd can be expressed as 

K d = kS. (44) 

AI2O3 smoke optical properties are used in evaluating the transfer equation since 
the smoke particles dominate the gas flow optical properties as discussed earlier in the 
report. Parry and Brewster [47] determined the optical properties of AI2O3 smoke 
produced by a solid propellant containing 20% aluminum by mass, burning at a 1 .8 MPa 

pressure. Optical thicknesses, Kq = 1 .28 at X = 632.8 nm and Kd = 1 .21 at X = 1064 
nm, were found for the 1 mm thick smoke region. Assuming a constant k in the 
transmission direction, Eq. 44 produces k = 1280 and k = 1210 for Xi and X2, 
respectively. 

Although the above k values were obtained from a solid propellant flame, they 
should provide a first approximation of the gas flow optical thickness in the model 
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combustion chamber. Given that k is a function of fume volume fraction, k actually may 
be greater for aluminum slurry combustion than for the above solid propellant due to the 
higher slurry aluminum mass percentage. As a conservative estimate, Parry and 

Brewsters' k values will be used in estimating a minimum gas flow optical thickness. 

In the model, chamber radius, R = 0.025 m, is the characteristic path length. 
Using this radius, the optical thickness is between 30 and 32. These values are much 
greater than Kq = 2.0, which is considered the lower limit for an optically-thick medium 
[51]. Therefore, as a first estimate, the gas flow can be treated as an optically-thick 
medium, in which local radiation heat transfer is only influenced by the immediate 
surroundings. Based on this optical thickness, a diffusion approximation [51] can be 
used to simplify the radiative transfer equation. The cylindrical Pi diffusion 
approximation [52] has been chosen since it provides greater accuracy than other 
approximations, yet is still simple to incorporate. A detailed development of this 
radiation model is presented in Ref. [29]. The final result of this analysis for the radiant 
heat transfer rate from the particles to the wall of a cylindrical chamber with radius, R, is 


q r (R) = 


4tt(I bf -I bw )ll(%) 


3 Uo % + 21,1/1; 



where 


e = 1 

V3(i-n 0 ) ’ 

and is the scattering albedo, 


(45) 


(46) 
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q s 

a + o s 


(47) 


where I 0 (R/£) and Ii (R/£) are the modified Bessel functions 



g c*r 

k =okir(K + v + 1) ’ 


where v equals 0 or 1 . 


(48) 


Since Q 0 , and therefore £, is a function of wavelength, X, the radiant heat 
transfer rate (Eq. 45) must be integrated over all wavelengths as follows: 


pi 
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q r (R) = 


(49) 
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If £(X) is assumed to be constant in a wavelength interval, Eq. 49 can be expressed as 
the following summation: 


q r (R)= X 

n=1 

where R i (T) is the fraction of total radiation emitted by a blackbody at 

A 'n+1~ A n v ' 

temperature T in the wavelength interval between X n +i and Xn. 

The absorption and scattering cross-sections, C a and C s are determined from 
Mie theory for each of the M wavelength regions in Eq. 49. The absorption and 

scattering coefficients, a and a s , are determined from the following relations [51]: 
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(50) 


a = C a N 

(51 a) 

<Ts = Cg N 

(51b) 


where N is the number density of fume particles, and it has been assumed that the 
AI 2 O 3 fume particles are uniform in size. 

Momentum Conservation - The gas-phase momentum equation in the 
combustion chamber is trivial, assuming a negligible chamber pressure gradient and no 
body forces. The negligible chamber pressure gradient condition should be accurate for 
current conditions, but does require the assumption of no wall frictional losses or 
pressure drops due to flow acceleration. However, the momentum equations governing 
the hydrocarbon, aluminum, and agglomerate AI 2 O 3 mass fluxes are significant. In a 
given slurry droplet, there is no slip between the hydrocarbon and the aluminum before 
hydrocarbon burnout, and no slip between the aluminum and the AI 2 O 3 agglomerate 
afterwards: therefore, the hydrocarbon aluminum and AI 2 O 3 agglomerate momentum 
equations in a given droplet size class can be reduced to a single momentum equation. 
Virtual mass and Bassett forces can be neglected since the particle density is much 
greater than the gas density. The particle momentum equation for a droplet size class 
can be expressed as 

F = ma, (52) 

where F is the drag force on a particle, m is the particle mass and a is the particle 
acceleration, dup/dt. Substituting for the drag force using a drag coefficient, Cp, and 
using the chain rule to relate dup/dt to dup/dx results in the following form of Eq. 52: 
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du p _ 3 Pg C D (u g -Up)|u g -u p 
dx 4 Ppd Up 

In this equation, p g is the gas density, u g the gas velocity, u p the droplet velocity, p p the 
droplet density, and d the droplet diameter. 

The drag coefficient is approximated as that of a sphere using the following 
correlation [52]: 



0.4, 


(54) 


where the Reynolds number, Re, is based on the slip velocity between the gas and the 
droplet. 

Slurry Combustion - The physics of slurry droplet combustion and a probable 
secondary atomization process are described above and in other references [18- 
20,24,37], After hydrocarbon burnout occurs in a droplet or the droplet fragments, the 
remaining agglomerate(s) of aluminum particles heats up, melts to form a single molten 
aluminum droplet, and burns. 

Cho and Takahashi's shell formation model [37] is used to predict the droplet 
diameter at which rigid shell formation occurs for each droplet size class in the 
combustor model. When running the combustor code, it is assumed that secondary 
atomization occurs when the droplet diameter reaches the predicted rigid-shell 
diameter. Although secondary atomization actually occurs some time after rigid shell 
formation, the time interval presently is not known and is therefore, neglected in the 
combustor model. However, future experimental efforts should provide an estimate of 
this time interval, and then the time interval will be included in the combustor model. 

Particle size distribution after secondary atomization is also currently an unknown 
to be determined from experimental measurements and is treated as a system variable. 
In the combustor model, a droplet undergoing secondary atomization is presently 
assumed to shatter into a specified number of equal-size secondary droplets. Defining 
the fragmentation ration, (5, as the number of secondary droplets produced per initial 
droplet, a new value of N/xj can be expressed as 

= P- (55) 

j,new z j.old 

Knowing the slurry mass flux of a given size class, m s iuny,j. a new size class droplet 
diameter, dj^ew. can be found from the mass conservation expression, 


N 

x 
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(56) 


m slurry,j “ 


N 


j,new 


*d 3 
6 )- new 


Because of the close proximity of the slurry droplets to each other, the droplets are not 
surrounded by individual flames. Therefore, droplet combustion can be modeled as an 
evaporation process. Detailed development of the hydrocarbon gasification problem are 
presented in Ref. [29]. The final mass balance for the droplet is 


dm droPi LH 

dt 


i^Viin c p.g^g T b) + •) 
C P.g hfv.LH+lj ^ 0 ' 1 c p,LH dT 


(57) 


which is used in the system mass and energy governing equations. 

As mentioned previously, aluminum combustion proceeds through two different 
oxidation mechanisms; the first is aluminum vapor oxidation, and the second is droplet 
surface condensation/oxidation. A simple aluminum combustion model, based on 
previous research [19], was developed and incorporated into the combustor code [29]. 
Because aluminum agglomerate heating was treated as a simple step change in 
agglomerate temperature from the hydrocarbon temperature to the aluminum boiling 
temperature [29], droplet heating was neglected in the aluminum combustion model 
[29]. Recently, modeling of the aluminum agglomerate heating process has been 
improved as described earlier in this report (cf. Eq. 39). However, agglomerate heating 
from the hydrocarbon temperature to the aluminum boiling temperature is now 
separated into two stages. In the first stage, the agglomerate is heated to a molten 
state at the aluminum melting temperature [Eq. 39]. In the second stage, droplet heat- 
up during combustion brings the droplet from the aluminum melting temperature to the 
aluminum boiling temperature. This droplet heat-up is represented in the aluminum 
combustion model as the heating of a thin droplet surface layer from the aluminum 
melting temperature to the aluminum boiling temperature. This representation closely 
parallels the surface layer heat-up in the hydrocarbon vaporization model [29], and is 
therefore, not detailed in this report. The final result is 


dnridrop,Ai 

dt 


47 tk g r s [n 
c p.g 


h fv,AI “ T l h fv,AI 2 03 


T b,AI 

+ j c p> A|dt 

T m A! 


+ 1 


(58) 


where q is the fraction of the aluminum oxide that condenses, or is formed, at the 
droplet surface. It should be noted that hfv > Ai 2 o 3 is not a true enthalpy of vaporization 

since AI2O3 does not exist in a gaseous state. Rather, it dissociates into 2AI + 3/2O2. 
Therefore, hf Vi Ai 2 o 3 is actually the enthalpy released by the chemical reaction, 2AI( g ) + 

3/2O2 -> Al2C>3(i), occurring at the droplet surface temperature. 
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Solution Method - Equations 32, 35, 53, 57, and 58 form the governing equation 
set and are numerically integrated in the axial direction using the IMSL integration 
routine DVERK [53] to determine rh g (rhh) g , u p , rti LH , m A |, and mAi 2 o 3 - These values 

are then used to determine the gas temperature, T g (x), the radiation heat flux, q(x), final 
AI 2 O 3 agglomerate diameter, and chamber burnout length. Initial conditions are 
supplied at the injector face. Reference [29] provides a description of the methods 
employed to determine the various required thermodynamic and transport properties. 


One-Dimensional Model Results 

The one-dimensional engine code was exercised using the chamber diameter, 
pressure, flow rates, and aluminum mass loading presented in Table I. These values 
were chosen to simulate Galecki's test conditions [43]. Presently, chamber length is 
varied to allow complete propellant combustion. In the present study, an arbitrary 
normalized droplet size distribution is used (Fig. 12). 

Table 1 

Model Operating Conditions 


Chamber Diameter 
Chamber Pressure 
Slurry Flow Rate 
Aluminum Loading 
Oxidizer Flow Rate 


0.0522 m 
690 kPa 
0.00142 kg/s 
60% 

0.00312 kg/s 


Secondary atomization effects on combustion chamber gas and droplet velocities 
can be seen in Figs. 13 and 14. The data in Fig. 13 were calculated assuming no 

secondary atomization while the data in Fig. 14 assume a fragmentation ratio, p, of 40. 

In both cases, smaller droplets equilibrate more rapidly with the gas velocity than large 
droplets, and the gas velocity increases along the combustor axis due to decreasing gas 
density as temperature rises and mass addition to the gas flow from the slurry droplets. 

As expected, droplet velocities are the same in both cases until secondary 
atomization occurs, after which, the slope of the Fig. 14 velocity profiles decreases 
sharply and the droplet velocities rapidly equilibrate with the gas velocity. This rapid 
equilibrium is caused by secondary atomization shattering the initial droplets into small 
secondary droplets which equilibrate with the gas velocity more quickly than the initial 
droplet could. 



NORM. 
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FIGURE 14. Gas and droplet velocities versus axial location for a 60 wt% Al 
slurry assuming no secondary atomization fragmentation ratio of 


Gas temperature and composition using a fragmentation ratio of five, are 
presented in Fig. 15. The jagged shapes of both the temperature and composition 
profiles arise from using only four droplet size classes instead of a continuous droplet 
size distribution. The inflection in the temperature profile between 0.05 and 0.15 m is 
caused by the large enthalpy transfer from the gas flow to heat the aluminum after 
hydrocarbon burnout. The composition plots show that the flow is fuel lean, as 
evidenced by excess O2 present after propellant burnout, and that CO2 and H2O 
undergo dissociation due to the high gas temperature. 

In Fig. 16, chamber burnout distance is plotted versus fragmentation ration, |3, to 
illustrate the benefits of secondary atomization. It is readily apparent that only slight 
secondary atomization is required to yield significant decreases in droplet lifetimes. 
Higher secondary atomization intensities, represented by larger fragmentation ratios, 
have a lesser effect on burnout length because droplet lifetime is inversely proportional 

to droplet surface area, which increases as (p) 273 . 

Final oxide agglomerate diameter as a function of secondary atomization 
intensity, p, is shown in Fig. 17. Similar to the trend seen in Fig. 16, small secondary 
atomization intensities significantly reduce final AI2O3 agglomerate diameters, with 
greater atomization intensities providing decreasing marginal reductions in diameter. 

However, final particle diameter is proportional to (1/p) 1/3 and not (1/p) 273 as was droplet 
lifetime. 

The above results demonstrate that the combustor model produces realistic 
results, given current model assumptions, and that only low secondary atomization 
intensities are required to significantly decrease chamber burnout length and final 
particle size. 


SUMMARY AND CONCLUSIONS 

Experimental and analytical methods were applied to obtain an improved 
understanding of the ignition and secondary atomization processes associated with 
aluminized hydrocarbon propellants. A numerical model of a rocket combustion 
chamber was also developed to help ascertain the influence of secondary atomization 
on combustor performance. The principal accomplishments and conclusions of our 
investigations are the following: 

1 . Apparatus and experimental techniques were developed to study the 
secondary atomization of slurry droplets in the 10-200 pm diameter size 
range. Specifically, non-intrusive optical methods were implemented to 
measure droplet size and velocity, and to detect aluminum combustion of 
free-flying droplets in a high-temperature environment. 

2. Models of slurry droplet disruption (secondary atomization) were applied to 
aluminum/liquid hydrocarbon propellants and show the following: 
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GAS TEMPERATURE (K) 



FIGURE 15. Mass flux of major gas species and gas temperature versus axial 
location. Data is for a 60 wt% Al slurry assuming a secondary 
atomization fragmentation ratio of 5. 
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FIGURE 17. 


Final oxide agglomerate diameter as a function of secondary 
atomization fragmentation ratio. Data is for a 60 wt% Al slurry. 





• Increasing solids mass fractions and decreasing ultimate particle diameter 
are predicted to significantly decrease the minimum slurry droplet diameter 
required to form a rigid shell, thereby reducing the minimum diameter 
required for secondary atomization. 

• Slurry droplets close to the minimum diameter for rigid shell formation 
should contain little or no liquid slurry after shell sealing and larger 
droplets should contain a significant amount of liquid slurry upon shell 
sealing. Therefore, secondary atomization may be less effective in the 
smaller droplets since only shell fragments, rather than secondary slurry 
droplets, are produced. 

• For small droplets (d 0 < 50 pm), shell stress is predicted to increase at 
essentially the same rate, regardless of the shell sealing time, as a 
consequence of constant shell thickness. In contrast, the stress in larger 
droplets was found to increase more rapidly for earlier sealing times since 
the shells are thinner at these shorter times, indicating that short sealing 
times may lead to rapid secondary atomization. 

Theoretical calculations, based on a critical rigid shell thickness of three 
ultimate particles, predict a larger minimum slurry droplet diameter (do.iim = 
34.7 pm) required for disruption than was observed in the experimental data 

(do.lim = 20-25 pm). This could be due to experimental uncertainties 
associated with the unknown slurry particle index of refraction or the critical 
shell thickness required for rigid shell formation. 

Calculated agglomerate ignition times were found to be in reasonable 
agreement with experimental results. 

A one-dimensional model of a rocket combustion chamber incorporating 
multiple droplet size classes, gas-phase chemical equilibrium, two-phase flow, 
hydrocarbon vaporization, aluminum combustion, simple secondary 
atomization, and radiation heat transfer was developed. Exercise of the 
model predicts the following: 

• Only moderate secondary atomization is required to effectively reduce 
overall propellant burn times and final AI 2 O 3 residual size. Greater 
secondary atomization intensities provide little additional benefit. 

• Preliminary results indicate that radiation losses to the chamber walls are 
less than 3% and should be even less in a full-size engine because of 
increases in combustion chamber diameter, and consequently, the flow 
optical thickness. 
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